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ABSTRACT. A growing amount of observational and theoretical evidence 
suggests that most main sequence stars are surrounded by disks of 
Gometary material. In this paper we investigate the dust production by 
comets in such disks when the central stars evolve up the red giant and 
asymptotic giant branch (AGB). Once released, the dust will be ablated 
and accelerated by the gas outflow and the fragments will become the 
seeds necessary for condensation of the gas. The origin of the 
requisite seeds has presented a well known problem for classical 
nucleation theory. This model is consistent with the dust production 
observed in M giants and supergiants (which have increasing 
luminosities) and the fact that earlier supergiants and most WR stars 
(whose luminosities are unchanging) do not have significant dust clouds 
even though they have significant stellar winds. Another consequence of 
the model is that the spatial distribution of the dust will not in 
general coincide with that of the gas outflow, in contrast to the 
conventional condensation model. A further prediction is that the 
condensation radius is greater than that predicted by conventional 
theory, in agreement with IR interf erometry measurements of a-Ori . 

1. INTRODUCTION 

The sources of the small dust grains observed by IRAS around Vega, 
g-Pic and numerous other old main sequence stars (Aumann, 1985; Backman 
and Gillett, 1987; Walker and Wolstencrof t , 1988) are believed to be 
disks of cometary material (Harper et al . , 1984; Weissman, 1984; Matese 
et al . , 1987). Modeling of the IR emission yields dust disk radii 
extending to several hundred AU (Gillett, 1986). The disk in g-Pic has 
been resolved optically and it extends to -1000 AU (Smith and Terrile, 

1987) . The mass of the small (^lOpm) grains directly observed is 

> O.O1M0 (e.g. Gillett, 1986), However, the lifetime of these small 
grains is less than the probable ages of the stars (Whitmire et al . , 

1988) and therefore more massive unobserved sources must also be 
present. Estimates of the total disk masses are highly uncertain since 
it must be assumed that a standard mass distribution index holds over a 
range from the observed small grains up to an assumed maximum mass. If 
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the maximum mass is taken to be comet size, total mass estimates range 
from ~ 1 O-3OOM0(Harpen et al . , 1984; Weisman, 1984; Gillett, 1986). 
Taking into account observational bias, it is likely that most main 
sequence stars (including the sun) possess disks of cool orbiting grains 
(Aumann, 1985; Backman and Gillett, 1987), implying that the phenomenon 
is not indicative of young systems as originally conjectured. Further, 
several F-type stars with 60pm excesses have age estimates of =2 X 10^ 
yr (Backman and Gillett, 1987). 

These observations are consistent with the theoretical expectation 
(Kuiper, 1951; Cameron, 1962) that a primordial disk of residual 
unaccreted planetesimals should exist Just beyond the planetary region 
of the Solar System. This disk has previously been invoked as the 
source of the observed steady state short period cornets (Fernandez, 

1980; Matese and Whitmire, 1986a, b) whose flux (Fernandez, 1980) and 
inclinations (Duncan et al., 1988) are inconsistent with their origin as 
captured long period comets. In addition to this residual flattened 
planetesimal disk there is also theoretical evidence suggesting the 
existence of an Isotropic massive (~100 ) inner Oort cloud extending 

inward to just outside the planetary region (Hills, 1981; Weissman, 

1985; Bailey, 1983). However, since the dust clouds around 
g-Pic and a-Psa are known to be flattened we assume that it is the 
residual planetesimal disk that dominates within several hundred AU of 
the star. 

In this paper we investigate the dust production by disk comets 
around intermediate mass stars (=1-10 Mq ) as they evolve into red giants 
and especially red supergiants (AGB phase). During the AGB phase the 
stellar luminosity can increase to IO^-IO^Lq , depending on mass. We 
will argue that if a modest population of comets (~1M^) exi stSgWi thin 
100 AU of these stars they can readily supply the seeds (> 10 A) 
necessary for dust condensation in the gas outflow, thereby mitigating 
the classical nucleation problem of statistically producing sufficient 
critical clusters in this environment (Donn and Nuth, 1985; Draine, 
1985). 

II. ANALYSIS 

Iben and Renzini (1983) have modeled the evolution of intermediate 
mass stars in the AGB phase. The time dependence of the stellar 
luminosity can be written as 

L(t) = L(0) exp(t/T ) (1 ) 

Li 

where 

L(0) = 26.5(M./M (2) 

and 

T, = 6.5 X 10^(M./M )“^*^^yr. (3) 

Li 1 
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In these formulae t = 0 at the beginning of the AGB phase, and is the 
initial zero-age main sequence stellar mass. These formulae are 
applicable for masses in the range IM^ ^ M. ^ IOM 0 , however for masses 
< should be set = 2 Mq . The stellar luminosity, radial distance 

R and black body temperature T are related by 

L(t) = (R/AU)^(T(t)/280K)^LQ (4) 

while the time-dependent black body temperature can be written 

T(t) = T(0) exp(t/x ) where the temperature time scale 

T' T L 

Since, in a standard cometary mass distri bution , most of the mass 
resides ig^the largest bodies, we consider comets of a single mass 
M = 10 g. The total mass loss rate for observed comets can be fit 
to the form (Ney, 1982) 

2/R 

dM/dt - -m ^exp(-1552K/T) (5) 


where the Boltzmann factor is appropriate for the activation energy of 
water released from clathrates and the coefficient k depends on specific 
surface absorptivity and activity, being smaller for older comets like 
Encke, This semiempirical formula is applicable for black body 
temperatures T up to =700K , which encompasses the temperatures of 
relevance in the present analysis. 


Equation (5) can be integrated to give 




1552K, _ ^ ,15b2K, 
T ^t(O) ^ 


640 yr _ ^.1/3^ 


( 6 ) 


where E^ is the Error function, f h m/M = fraction of original mass 
surviving. Th^ nynjerical factor was evaluated by taking 
k = 3*7 X 10 g s (comet Encke). We shall assume that the second 
term is negligible compared to the first, an adequate approximation 
after several luminosity timescales have elapsed. Eq. (6) gives the 
value of T when M has been reduced to 

max 

We now locate ^that position in the disk where sublimation is 
maximal. Setting M = 0 we find 


f1/3 

exp (-1552K/T ) 640 yr p 

= ( 7 ) 

1552K/T Xt 2 

P b 

where Tp is the black body temperature where the comet is undergoing 
peak sublimation and fp is the surviving mass fraction at peak 
sublimation. T is also a rough lower limit to the black body 
temperature at the inner edge where f = 0. An adequate approximate 
solution to this equation is found to be 
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l860K/ln(T /400 yr). 

i-j 


( 8 ) 


T 


Having obtained the temperature at which comets are undergoing peak 
sublimation rates we next determine the disk location where this 
occurs and take this as an adequate (upper bound) approximation to the 
inner edge, Rp is found fran 

Rp(t) = (L(t)/Lg )V(Tp/280K)^AU (9) 

or 

Rp(t)/R^(t) = *4(T^/Tp)^ (10) 

where * refers to stellar parameters. In Table 1 Tp, fp and Rp/R* 
tabulated for several stellar masses. Rp/R* is given for two stellar 
temperatures corresponding to a -Ori (oxygen supergiant) and IRC 10216 
(carbon supergiant). 


Table 1 


Rp/R* 


Mi /Mg 

Tp(K) 

fp 

T^e=3300K 

Tjt=2300K 

2 

1 94 

0.33 

145 

70 

4 

270 

0.34 

75 

36 

6 

341 

0.35 

47 

23 

8 

423 

0.36 

30 

15 

1 0 

509 

0.37 

21 

1 0 


Finally, if one approximates that sublimation occurs only at the 
location of the peak sublimation 

(t)- (2Tj_)-'(BdM/dR)|j 

P P 

For any standard radial power law RdM/dR ^ AM(i^ -► 2R) in order of 
magnitude. Therefore, if there is^--1M^of comets available we should be 
releasing'it at a rate 

/yr for M. = 2Mg, 

dM ^/dt ~ AM/t,^ <■" (12) 

net L > 

' /yr for M. = lOM^ 

To order of magnitude this result is of general applicability and can 
also be used, for example, to estimate the dust production rate around 
red giants as well as red supergiants. 
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III. DISCUSSION 


Although these dust production rates are incapable by themselves of 
explaining the observed dust clouds aromd most M giants and supergiants 
they can readily supply the seeds (> 10A) necessary for dust 
condensation in the gas outflows. Seeds of radii ^lOA whicn 
subsequently grow to ~ 0.1 pm represent a dust mass enhancement factor 
of ^10^. The required seeds should be readily produced by ablation in 
the dense gas outflow, sputtering, grain-grain collisions, or they could 
simply correspond to the lower end of a standard mass distribution of 
the sublimating particles. The seeds will simultaneously be accelerated 
until there is little relative velocity between the seeds and the gas 
outflow, at which point condensation may occur rapidly. We note that 
the tefnperature T„ in Table 1 is the black body tamperature at the 
location of those^comets which are undergoing ^peak sublimation, but the 
temperature of the small condensing grains (S: 0.1pm) that are actually 
observed can be significantly higher (e.g. 1000K). 

The inner radius of the best studied oxygen supergiant a Ori has 
been resolved with IR interf erometry and found to be = 35 - 60 R» 
(Bloemhof, 1984; Howell et al., 1981), consistent with the proposed 
model but much larger than expected by conventional condensation theory 
(Draine 1985). In other cases it is often assumed that T = 1000K which 
implies’an inner edge =5R* for black body radiators (Tielens, 1988). 

The proposed model is consistent with the fact that M giants and 
supergiants have significant dust clouds while earlier supergiants and 
most WR stars do not, even though they have significant mass loss via 
stellar winds. It is only the cool giants and supergiants that have 
rapidly Increasing luminosities , .and therefore short time scales , as 
they evolve. Most of the comets within a few hundred AU should be 
concentrated in a disk. Thus, we expect the sublimated seeds and 
condensed dust to be preferentially concentrated near the disk. Dust 
disks are compatible with the observed asymmetries in the IR emission in 
IRC 10216 and other cool giants (e.g. Zuckerman, 1980; Beckwith, 

1985). A prediction of the model is that the spatial distribution of 
the dust will not in general coincide with that of the gas outflow, in 
contrast to conventional condensation theory. Although R is uncertain 
by a factor of = 2, the model predicts that the inner disk edges will 
tend to be further out than suggested by conventional condensation 
theory, especially for the lower mass stars. 
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